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Ri jke 管 热 声 不 稳定 非 线 性 动力 学 研究 
-分 盆 分 析 


zy, HOC*, A Ase 
(西北 工业 大 学 ， 燃 烧 、 流 动 和 热 结构 国家 级 重点 实验 室 ， 西 安 710072) 
摘 要 : 建立 了 水 平 Rijke 管 热 声 模型 ， 利 用 Galerkin 方法 对 控制 方程 组 进行 展开 ， 实 现 数值 求解 。 确 定 了 
Galerkin 模 态 收敛 阶 数 为 10， 利 用 非 线性 动力 学 理论 对 系统 进行 分 岔 分 析 。Rijke 管 热 声 系 统 分 命 行为 属于 亚 临 界 
Hopf 分 岔 。 系 统 稳定 性 区 域 分 为 全 局 稳定 、 全 局 不 稳定 及 双 稳 态 区 域 。 获 得 了 无 量 纲 加 热 功率 开 ， 加 热 器 位 置 闵 ， 
阻尼 系数 cl 和 时 间 延 迟 fz 等 参数 的 分 命 图 谱 ， 发 现 加 热 器 位 置 x 的 分 命 图 谱 存 在 两 个 Hopf 分 岔 点 。 在 线性 不 稳定 
区 域内 ， 振 荡 幅 值 随 时 间 延 迟 t 的 增 大 呈现 先 增 大 后 减 小 的 趋势 。 
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Rijke Tube Thermoacoustic Instability Nonlinear 
Dynamics Study—Bifurcation Analysis 


FENG Jian-Chang AO Wen LIU Pei-Jin 
(Science and Technology on Combustion, Internal Flow and Thermal-structure Laboratory, 


Northwestern Polytechnical University, Xi’an 710072, China) 


Abstract: The thermoacoustic model of a horizontal Rijke tube is established. The governing equations 
are expanded and solved by using the Galerkin method. The convergence order of Galerkin modes is 
determined as 10. The bifurcation analysis of the system is carried out by using nonlinear dynamics theory. 
The bifurcation behaviors of the horizontal Rijke tube thermoacoustic system is found to be subcritical Hopf 
bifurcation. The system stability regions are divided into globally stable region, globally unstable region and 
bistable region. The bifurcation diagrams of the non-dimensional heater power (K), the position of the heater 
(xp), the damping coefficient (cı) and the time delay (t) are obtained, and the bifurcation diagram of the 
position of the heater (xy) has two Hopf bifurcation points. In linear unstable region, the amplitude of the 
oscillation firstly increases and then decreases with the increase of time lag (7). 
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燃烧 不 稳定 ‘combustion instability) 是 指 燃烧 
器 中 压力 和 速度 的 大 幅度 自 激 振荡 ， 这 一 振荡 带 来 
的 循环 负荷 会 造成 结构 过 度 振 动 和 燃烧 室 壁 面 的 热 
量 过 度 传递 等 问题 ， 对 结构 组 件 造成 致命 的 破坏 。 
长 久 以 来 ， 燃 烧 不 稳定 一 直 都 是 火箭 发 动机 ， 燃 气 
轮机 ， 航 空 发 动机 和 电站 锅炉 等 许多 高 性 能 燃烧 装 
置 面 临 的 重大 挑战 。 热 声 不 稳定 〈thermoacoustic 
instability) 是 燃烧 不 稳定 的 最 为 主要 的 一 种 形式 ， 
一 般 认 为 ， 热 声 不 稳定 是 由 于 不 稳定 放 热 和 声场 所 
荡 的 耦合 引起 的 ， 将 引发 接近 固有 声场 频率 的 大 
幅 压 力 振荡 ， 其 本 质 是 热源 的 放 热 与 声场 环境 存在 
正 反馈 机 理 口 。 
燃烧 不 稳定 现象 背后 所 蕴含 的 机 理 十 分 复杂 ， 
包含 燃烧 ， 流 动 和 声场 的 相互 耦合 作用 。 火 焰 面 的 
变化 ， 当 量 比 波动 ， 流 体 动 力学 不 稳定 导致 的 涡 脱 
落 以 及 燃料 的 雾 化 及 薰 发 过 程 中 的 振荡 都 可 能 驱动 
燃烧 不 稳定 。 而 燃烧 室 几 何 结构 ， 人 燃料 /空气 混合 
度 ， 燃 料 的 类 型 和 构成 ， 以 及 工作 环境 (注入 空气 
温度 ， 人 燃烧 室 压强 ， 当 量 比 ， 涡 流 强 度 ， 贫 油 预 混 
等 ) 都 将 对 燃烧 不 稳定 产生 影响 叫 。 
流动 和 火焰 的 非 线性 动力 学 特性 是 燃烧 不 稳定 
中 最 为 重要 的 现象 之 一 。 燃 烧 室 中 的 流动 和 火焰 特 
性 随 着 控制 参数 的 变化 会 发 生 剧 烈 改变 ， 尤 其 是 当 
参数 变化 经 过 分 岔 点 时 。 燃 烧 过 程 本 身 可 能 会 ， 抑 
或 不 会 出 现 分 岔 现 象 ， 但 一 旦 发 生 分 岔 ， 非 线性 行 
为 在 许多 燃烧 装置 中 将 占据 主导 地 位 B4。 系 统 参数 
的 微小 扰动 有 可 能 会 引起 分 岔 ， 导 致 系统 由 稳 态 
(小 振荡 或 无 振荡 ) 转换 为 不 稳定 状态 〈 大 振幅 的 
极限 环 )。 相 反 地 ， 在 相同 的 临界 参数 值 时 ， 系 统 
由 不 稳定 状态 转变 为 稳 态 则 不 一 定 会 发 生 ， 这 是 由 
于 滞后 现象 的 存在 。 
燃烧 不 稳定 中 出 现 的 分 贫 ， 汪 后 ， 极 限 环 和 触 
发 等 现象 均 可 以 用 动力 系统 理论 来 解释 和 描述 ， 这 
。 通 过 非 线 性 动力 学 
理论 ， 可 对 非 线 性 系统 的 特征 进行 深入 的 刻画 、 分 
析 和 诊断 。 现 有 研究 大 部 分 以 燃气 轮机 和 火箭 发 动 
机 出 现 的 热 声 不 稳定 问题 为 背景 ， 研 究 对 象 则 通常 
为 各 种 简化 的 燃烧 器 或 实际 的 燃气 轮机 和 发 动机 。 
其 中 ，Rijke 管 是 研究 热 声 现象 最 方便 最 典型 的 实验 
系统 ， 也 是 研究 燃烧 室 火 焰 区 燃烧 与 声场 的 相互 作 
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用 导致 的 不 稳定 机 理 及 控制 的 基础 ， 因 此 有 许多 研 
究 是 针对 Rijke 管 来 开展 的 58。 不 稳定 放 热 作为 热 
声 不 稳定 的 激励 ， 是 不 稳定 发 生 的 一 个 关键 环节 。 
热 声 系统 中 的 能 量 输入 可 以 有 多 种 形式 ， 例 如 燃烧 
室 中 放 热 的 化 学 反应 ， 电 热 丝 ， 电 热 薄 膜 Cot 
gauze) %0], Jahnke 和 Culick19 最 早 将 现代 动力 系 
统 理论 引入 燃烧 不 稳定 研究 之 中 ， 以 分 析 非 线性 燃 
烧 不 稳定 。 他 们 提出 了 一 种 延 拓 算法 ， 可 对 各 稳 态 
和 极限 环 的 解 进行 体系 化 高 效 计算 ， 由 此 可 得 到 稳 
态 和 极限 环 下 的 稳定 性 ， 通 过 分 命 分 析 可 以 判定 不 
稳定 发 生 点 。 不 稳定 发 生 点 处 的 声波 的 性 质 取决 于 
分 岔 的 类 型 和 分 岔 点 处 极限 环 的 性 质 。 
Ananthkrishnan 等 [1 针对 燃烧 室 中 的 不 稳定 运动 建 
立 了 降 阶 模型 ， 并 求解 分 析 了 触发 和 极限 环 性 质 。 
分 析 表 明 ， 对 于 轴 向 振荡 模 态 ， 一 阶 不 稳定 需要 至 
少 四 阶 Galerkin 级 数 展 开 ， 而 二 阶 不 稳定 则 至 少 需 
要 八 阶级 数 展开 ， 才 能 确保 计算 的 准确 性 。 
Balasubramanian02 研 究 了 非 正 则 性 在 简单 热 声 系统 
中 的 作用 和 影响 ， 通 过 推导 水 平 Rijke 管 的 声场 控 
制 方程 并 求解 ， 重 点 讨论 非 正 则 性 和 非 线 性 效应 带 
来 的 结果 。Subramanian03] 分 析 了 热 声 系 统 中 亚 临 界 
分 岔 的 性 质 ， 以 及 随后 导致 的 双 稳 定性 问题 。 
Waugh 等 0415 研 究 了 热 声 系统 中 的 触发 和 流体 力学 
中 的 向 满 流 的 旁 路 转换 〈“ 请 越 型 ”转换 ) 的 比 
拟 ， 二 者 机 理 均 是 小 扰动 导致 系统 变 为 大 幅度 自持 
振荡 ， 尽 管 系统 是 线性 稳定 的 。 研 究 证 明了 某 些 类 
型 的 噪声 更 能 触发 振荡 。 

现 有 研究 对 Rijke 管 热 声 不 稳定 非 线 性 动力 学 
有 一 定 认 识 ， 但 对 系统 分 岔 性 质 的 认识 尚 不 够 清 
晰 ， 对 于 无 量 纲 加 热 功 率 、 加 热 器 位 置 、 阻 尼 系 数 
等 系统 参数 对 分 岔 的 影响 规律 缺乏 了 解 。 本 文 建立 
T Rijke 管 热 声 不 稳定 的 简化 模型 ， 得 到 其 一 维 声 
场 的 控制 方程 ， 并 对 其 进行 了 求解 ， 通 过 对 Rijke 
管 的 热 声 不 稳定 性 进行 分 岔 分 析 ， 加 深 对 Rijke 管 
热 声 不 稳定 产生 机 理 和 振荡 现象 的 认识 。 
1 热 声 不 稳定 模型 及 求解 
1.1 Rijke 管 模型 

本 文 所 研究 的 水 平 Rijke 管 模型 如 图 1 所 示 。 
在 该 模型 中 ， 管 长 工 是 固定 不 变 的 ， 加 热 器 的 位 置 
xy 的 值 可 以 通过 改变 加 热 器 在 管道 内 的 位 置 来 改 
变 ， 热 源 为 电热 丝 加 热 装置 ， 在 实验 条 件 下 其 加 热 


强度 可 以 通过 改变 通过 的 电流 来 改变 。 


Heat source 


图 1 KF Rijke 管 的 结构 示意 图 ， 来 流 平均 速度 为 4 ， 管 长 
为 L， 热 源 与 来 流 开口 端 距离 为 x。 


Fig.1 Structure schematic of a horizontal Rijke tube. The mean 


flow velocity is u , the length of the duct is Z and the distance 


between the heat source and the incoming flow open end is x» 


1.2 控制 方程 


忽略 平均 流量 与 温度 梯度 的 影响 ， 控 制 管内 声 
场 的 动量 和 能 量 方程 如 下 [7]:; 
_Ou Op’ 
Pa ag SD 
pet +op'=(y-NO'() (2) 
t Ox 


为 了 简化 和 方便 分 析 ， 对 式 〈1) 和 2) 进行 


无 量 纲 化 处 理 ， 得 到 无 量 纲 化 的 系统 控制 方程 组 : 
十 0 (3) 
Ot Ox 


P ym, rep = -DOOS 4) 
t Ox 


其 中 无 量 纲 化 的 尺度 如 下 : 


(5) 
式 (1) - (5) 中 ，x 是 沿 轴 向 方向 的 距离 ，t 是 时 
ll], p 为 管内 的 平均 压强 ，wu' 为 声波 速度 ，p' 为 
声波 压强 ， 2 为 流体 密度 ，y 是 介质 的 比 热 比 ， 声 
HN c, tM, 为 平均 流 马赫 数 。 此 外 ，< 上 为 阻尼 系 
BH, Ot) 为 单位 面积 上 由 于 电 加 热 器 产生 的 热 释 放 
速率 波动 ， 5 为 标准 狄 拉 克 分 布 ，xy 为 热源 位 置 。 
符号 “” 和 ““ ”分别 表示 有 量 纲 量 和 平均 量 。 


阻尼 可 以 用 下 式 表示 ; 
al I wi 
Dera a (6) 


参数 c Me, 保持 常数 不 变 ，wi = jz 为 第 j 阶 管内 声 


a eT W 
CM aX | Vo 


模 态 的 无 量 纲 频率 。 
热源 模型 采用 Heckl 经 验 公 式 09: 
ó) = 
DEAT. =T), 
SV3cop ACP ents owed g 


KH, L, oe ee T, 
为 电热 丝 的 温度 ， 了 代表 周围 空气 的 平均 温度 ，8 
为 管道 截面 积 ，4 代表 空气 的 热传导 率 ，c, 表示 恒 
定 体积 内 单位 质量 空气 的 比热容 ，2 为 气体 的 平均 
密度 。 由 于 热 惯性 的 存在 ， 传 热 和 流动 速度 之 间 有 
一 个 时 间 延 迟 z 。 需 要 注意 的 是 ， 此 热源 模型 只 适 
用 于 “紧凑 型 ”热源 的 描述 ， 即 热源 厚度 相对 波长 
而 言 可 以 忽略 不 计 。 
将 式 (7) FRA FES 


量 方程 《4) 中， 则 能 量 方程 


‘al 


可 改写 为 
Op' ôu' 
Bt YM + Ep! =l- x,y -)) 
2L,(T,-T) _d,-| fll; E 
x sie, TAC, P 5 i T) J 


(8) 
1.3 Galerkin 展开 
求解 Rijke 管 装置 热 声 问题 的 常规 方法 是 使 用 
Galerkin 投影 和 模 态 扩展 。Galerkin 方法 的 原理 是 
任何 域 中 的 函数 都 可 以 表达 成 该 域内 展开 函数 的 登 
加 。 基 础 函数 的 选取 原则 是 它们 须 满足 边界 条 件 ， 
且 基 础 函数 的 选择 不 是 唯一 的 。 在 这 样 的 方法 中 ， 
从 满足 边界 条 件 的 空间 基 函 数 来 看 ， 空 间 和 时 间 上 
变化 的 压力 和 速度 信号 被 扩展 ， 这 些 基 函 数 的 选择 
不 是 唯一 的 。 本 文选 择 的 基础 函数 是 任意 的 ， 并 非 
系统 的 特征 函数 ， 而 是 线性 化 后 的 系统 自 共 力 
(self-adjoint) 部 分 的 特征 函数 。 因 此 ， 速 度 场 和 
压力 场 可 以 表示 成 管道 固有 模 态 的 形式 ; 


N 
= >》 Ui cos(k x) (9) 
j=l 
N 
p'=yM >" P, sin(k;x) (10) 
j=l 


N 
4 


P, ke, = jz 为 第 j 个 管道 模 态 的 无 量 纲 波 


长 。 在 极限 No 条件 下 ， 这 些 基 函数 形成 一 个 完 
整 的 基础 系 ， 显 然 ， 本 文选 取 的 展开 函数 满足 边界 
条 件 。 
将 式 (9) 和 C10) 替换 到 式 G) Fil (8) 

中 ， 并 投影 (projecting along) 到 基础 函数 系 ， 得 
到 : 


U,+k,P,=0 (11) 


一 3 sin(k ;xy) 


(12) 


P, +26;@;P,-k;U; 


A 


其 中 ， 无 量 纲 加 热 功率 为 : 


4( -DZ _ 
K-02 Dh T-P hiep a (13) 
yMaSc, pv3 2 


上 述 完 全 耦合 的 Galerkin 模型 可 以 通过 线性 化 
延迟 进一步 简化 ， 该 方法 只 适用 于 满足 z<7 的 
Galerkin 模 态 下 ， 其 中 7 =2/j 为 第 j 阶 Galerkin 模 
态 的 时 间 周 期 。 在 满足 这 种 假设 的 条 件 下 ， 延 迟 项 
可 以 写成 : 


1 N 
一 十 》 cos(k.x,)U,(t-—t 
2, (kx, U,(t-7) 


U,(t—7) xU, + th, P (14) 
通过 线性 化 延迟 ， 得 到 最 终 Galerkin 展开 的 常 
微分 方程 组 : 


U,+k,P,=0 GED 


= [Fors 


(j =1,2,....N.) (16) 


P,+26;0;P; kU; 


A 


其 中 ，k; = jz 为 第 j 阶 管内 声 模 态 的 无 量 纲 波 数 。 
1.4 Galerkin 模 态 收敛 

在 极限 V > REF, A O) 和 (10) 这 些 
基 函 数 形成 一 个 完整 的 基础 系 ， 但 是 考虑 到 计算 的 
可 行 性 ， 我 们 只 能 使 用 有 限 数 量 的 模 态 数目 。 精 确 
捕捉 到 系统 的 线性 和 非 线 性 行为 所 需 的 Galerkin 模 
态 的 数目 被 称 为 模 态 收敛 。 利 用 Matlab 编程 实现 了 
对 式 (15) 和 (16) 的 数值 求 ， 对 不 同 数目 的 
Galerkin 模 态 下 的 系统 演化 进行 比较 。 对 于 在 线性 


N 
3+ PU +7k,P) 
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不 稳定 区 域 ， 具 有 不 同 数量 声学 模 态 的 该 系统 的 时 
间 演 化 的 比较 见 图 2。 可 以 得 到 1 阶 Galerkin 声学 
模 态 下 极限 环 的 振幅 和 相位 与 10 阶 Galerkin 声学 
模 态 下 极限 环 的 振幅 和 相位 差别 较 大 ， 但 随 着 
Galerkin 模 态 数目 的 增加 ， 极 限 环 的 振幅 和 相位 差 
别 逐 渐 减 小 。 根 据 图 2， 我 们 可 以 看 出 ， 当 取 9 阶 
和 10 阶 声学 模 态 得 到 数值 解 的 差异 是 几乎 可 以 忽 
略 不 计 ， 也 再 次 说 明了 10 阶 声学 模 态 Galerkin 投 
影 的 收敛 性 。 对 于 10 阶 声学 模 态 ， 增 加 一 阶 
Galerkin 模 态 所 改变 的 声波 速度 极限 环 振幅 小 于 
1.4%021。 因 此 ， 在 本 文 以 后 的 所 有 计算 中 ， 我 们 采 
用 10 Bit Galerkin 模 态 的 模型 来 进行 线性 和 非 线性 
稳定 性 分 析 以 确保 收敛 性 。 
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2 不 同 声学 Galerkin 模 态 数目 下 系统 随时 间 演 化 的 比较 ， 
系统 参数 为 : K=1, cı=0.1, c=0.06, x=0.3;, 1=0.2. 


Fig.2 Comparison of time evolution of the system with different 


acoustic Galerkin mode numbers. The parameter values of the 


system are K=1, c\=0.1, co=0.06, x/=0.3 and t= 0.2. 
2 线性 稳定 性 分 析 


为 研究 无 穷 小 扰动 对 系统 稳定 性 的 影响 ， 首 先 
对 系统 进行 线性 稳定 性 分 析 。 如 果 系 统 逐 渐 远 离 稳 
定 状态 ， 则 系统 最 终 是 不 稳定 的 ， 如 果 系 统 逐 渐 接 
近 稳 定 状 态 ， 那 么 系统 最 终 达 到 稳定 状态 。 这 是 指 
局 部 分 析 系 统 在 稳定 状态 附近 时 的 稳定 性 变化 。 男 
一 方面 ， 非 线性 稳定 性 分 析 是 研究 有 限 振幅 扰动 对 
系统 的 影响 效果 ， 并 用 于 表征 所 得 到 的 系统 渐 近 状 
态 。 由 于 系统 平衡 解 失 稳 ， 系 统 的 行为 发 生 改 变 的 
临界 点 称 为 分 岔 点 。 一 旦 一 个 分 岔 点 根据 某 个 参数 
被 确定 ， 该 分 岔 点 本 身 继续 与 系统 另 一 个 相关 参数 
变化 。 分 岔 点 分 支 提供 了 线性 稳定 的 边界 ， 该 边界 
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可 以 分 离 相关 参数 空间 中 线性 稳定 和 线性 不 稳定 的 
区 域 。 这 个 稳定 边界 是 在 所 有 自由 参数 变化 空间 的 
个 超 曲 面 ， 但 比较 方便 的 是 我 们 可 以 表示 几 个 适 
当 的 二 维 投影 曲线 。 
在 Rijke 管 中 ， 系 统 参数 有 加 热 器 的 无 量 纲 加 
热 功 率 、 加 热 器 的 位 置 、 阻 尼 系 数 以 及 时 间 延 迟 ， 
但 是 在 实际 实验 操作 条 件 下 ， 可 以 精确 改变 的 参数 
只 有 无 量 纲 加 热 功率 和 加 热 器 的 具体 位 置 ， 时 间 延 
迟 和 阻尼 系数 都 是 根据 其 他 条 件 来 改变 的 。 我 们 以 
加 热 器 位 置 和 时 间 延 迟 项 的 线性 稳定 性 边界 为 例 ， 
对 加 热 器 位 置 和 时 间 延 迟 变 化 时 的 线性 稳定 性 区 域 
和 线性 不 稳定 区 域 进 行 分 析 。 一 个 典型 的 加 热 器 的 
位 置 和 时 间 延 迟 的 线性 稳定 边界 变化 如 图 3 所 示 ， 
显示 了 对 于 确定 的 阻尼 和 加 热 器 功率 等 固定 系统 参 
数值 ， 对 于 加 热 器 位 置 (xy) 的 一 个 选择 范围 ， 系 
统 是 否 为 线性 不 稳定 则 取 诀 于 时 间 延 迟 Cr), RZ 
亦 然 。 对 于 较 小 或 合理 大 的 r ， 比 如 rz<0.15 和 
T>0.85 ， 加 热 器 在 任何 位 置 时 ，Rijke 管 系统 都 是 
线性 稳定 的 ， 即 对 于 小 扰动 ， 系 统 都 最 终 趋 于 稳定 
状态 。 只 有 在 0.15 <rz < 0.85 范围 内 ， 系 统 平衡 解 的 
稳定 性 取决 于 加 热 器 的 位 置 。 


jm 


f Linearly stable region 
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图 3 加 热 器 位 置 六 和 时 间 延 迟 r 的 线性 稳定 性 边界 , 系统 其 
也 参数 为 : K=0.8, cı=0.1, c2=0.06. 


Fig.3 Linear stability boundary between x, and t. The other 


parameter values of the system are K=0.8, cı=0.1 and c2=0.06. 
3 分 叉 分 析 

分 命 分 析 是 研究 分 命 现象 的 特性 和 产生 机 理 。 
分 兮 现象 是 指 对 于 某 些 完全 确定 的 非 线 性 系统 ， 当 
系统 的 某 一 参数 4w 连续 变化 到 某 个 临界 值 IN, A 
统 的 全 局 性 性 态 〈 定 性 性 质 、 拓 扑 性 质 等 ) SRE 
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突然 变化 。 u, 称 为 参数 u WARTS. RA 
稳定 状态 失去 稳定 性 出 现 孤 立 的 极限 环 周 期 解 的 分 
BRN Hopf 分 岔 。 极 限 环 的 新 兴 分 支 的 稳定 性 决 
定 了 Hopf 分 支 的 性 质 或 类 型 。 这 两 种 Hopf TAK 
型 图 4 所 示 。 如 果 极 限 环 是 不 稳定 ， 如 图 4 (a) 所 
示 ， 不 稳定 极限 环 分 支 形成 一 个 双 稳 态 区 域 ， 在 该 
区 域内 稳 态 解 对 于 小 扰动 是 稳定 的 而 对 于 大 扰动 是 
不 稳定 的 。 然 而 ， 不 稳定 极限 环 的 这 一 分 支 可 能 会 
发 生 折 县 分 贫 和 稳定 平衡 。 对 于 小 于 折 肢 点 的 参数 
值 ， 稳 定 解 对 于 任何 幅度 的 扰动 都 是 稳定 的 ， 因 此 
系统 是 全 局 稳定 的 。 因 此 在 亚 临 界 分 岔 情 况 下 ， 线 
性 (局 部 ) 和 非 线 性 (全 局 ) 的 稳定 性 边界 是 不 同 
的 。 这 种 分 侈 行为 被 称 为 亚 临 界 Hopf 分 岔 。 但 如 
果 极 限 环 是 稳定 的 如 图 4 b) 所 示 ， 系 统 平稳 地 从 
稳定 的 稳 态 解 发 展 为 一 种 具有 递增 极限 环 振幅 不 稳 
定 的 稳 态 解 。 这 种 分 岔 类 型 称 为 超 临界 Hopf 分 
a> 


iO o 


(b) 


Measure 


Parameter 
图 4 (a) 亚 临 界 和 Cb) 超 临 界 Hopf 分 岔 附近 参数 变化 时 
对 应 的 分 命 行为 031。 


Fig.4 Bifurcation behavior of a measure for the variation of a 


parameter near (a) Subcritical and (b) Super-critical Hopf 


bifurcation. 

可 通过 定性 分 析 方 法 判定 亚 临 界 和 超 临 界 Hopf 
DE. (LE X)* 形式 的 源 项 非 线性 分 贫 的 性 质 可 以 
通过 把 非 线性 项 并 展开 成 级 数 ， 并 放弃 高 阶 项 来 判 
Wie EIR (14. X) 表达 式 的 二 项 扩展 的 结果 为 下 二 


a(a-I(a-2) ps 
3! 


- ssi 
(17) 
在 上 述 表 达 式 中 ， 一 阶 项 
aX 和 三 阶 项 w(w-1(w-2) /3! 的 符号 是 相同 
Pr 
aļ(a-1)(a-2) X /3! 这 两 项 的 符号 是 不 同 的 。 一 阶 
re eh 每 当 这 
两 项 具有 相同 的 符号 时 ， 分 命 是 亚 临 界 分 命 类 型 
的 ， 而 当 这 两 项 具有 不 同 的 符号 时 ， 分 岔 为 超 临界 
分 岔 类 型 的 。 对 于 Rijke 管 模型 中 的 热 释放 速率 波 
动 ，w=J/2 ， 这 意味 着 该 模型 将 出 现 亚 临界 类 型 
Hopf 7} 7 © 
在 Rijke 管 热 声 系统 中 ， 可 变 的 系统 参数 有 有 
加 热 器 的 无 量 纲 加 热 功 率 久 、 加 热 器 的 位 置 w%、 阻 
尼 系 数 c 以 及 时 间 延 迟 r 。 分 别 对 系统 内 无 量 纲 加 
ADK MARKE xr 阻尼 系数 c 以 及 时 间 
延迟 z 的 变化 进行 分 岔 分 析 。 
3.1 加 热 器 功率 的 影响 
在 Rijke 管 热 声 系统 中 的 加 热 器 为 电热 丝 加 
热 ， 通 过 改变 K 值 来 探究 加 热 器 功率 对 系统 产生 的 


影响 。 无 量 纲 加 热 器 功率 可 以 通过 增加 提供 给 加 热 
器 的 功率 来 达到 ， 其 增 大 代表 着 系统 驱动 力 的 增 
加 ， 增 加 的 驱动 力 使 得 系统 变 得 更 加 不 稳定 。 因 
此 ， 对 于 较 小 的 玉 值 ， 系 统 的 平衡 是 稳定 的 ， 在 所 
有 初始 扰动 下 ， 都 会 衰减 渐 近 接近 于 零 。 增 加 天 值 
将 降低 流动 的 稳定 性 裕 度 ， 在 分 岔 点 附近 ， 系 统 演 
变 为 线性 不 稳定 ， 从 而 导致 管内 的 振荡 流 模式 。 图 
5 为 计算 得 到 的 无 量 纲 加 热 器 功率 〈 天 ) 的 变化 对 
系统 演化 的 影响 。 空 心 圆 表明 不 稳定 解 ， 实 心 圆 代 
表 稳定 解 。 由 图 可 知 靠近 Hopf 分 兮 点 的 小 振幅 极 

限 环 是 不 稳定 的 ， 这 些 不 稳定 的 极限 环 与 稳定 平衡 
解 共存 。 这 些 极限 环 的 不 稳定 分 支 进一步 经 历 一 个 
折 礁 点 分 岔 获得 稳定 。 据 此 可 判定 分 岔 属于 亚 临界 
类 型 ， 证 实 了 之 前 的 结论 。 在 其 他 系统 参数 一 定 的 
条 件 下 ， 系 统 随 时 间 的 演化 行为 取决 于 天 值 的 变 

化 。 在 玉 <0.515 时， 系统 对 于 任何 幅度 的 扰动 ， 都 
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最 终 趋 于 零 ， 此 区 域 为 全 局 稳定 区 域 。 随 着 天 值 的 
不 断 增 大 ， 当 0.515 < 天 < 0.615 时 ， 系 统 有 三 种 可 
能 的 状态 ， 稳 定 状态 ， 不 稳定 极限 环 以 及 稳定 极限 
环 。 根 据 初 始 条 件 的 不 同 ， 系 统 可 能 最 终 进入 稳定 
状态 或 者 周期 振荡 ， 该 区 域 被 称 为 双 稳 态 区 域 。 当 
K >0.615 时 ， 系 统 进入 全 局 不 稳定 区 域 ， 对 于 任意 
小 幅度 的 初始 扰动 ， 系 统 最 终 都 会 达到 极限 环 状 
态 。 通 过 无 量 纲 加 热 功 率 的 分 爹 图 也 可 以 证 明 无 量 
纲 加 热 功 率 的 增 大 是 系统 产生 热 声 不 稳定 的 主要 原 
因 。 
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Fig.5 Bifurcation plot for variation of non-dimensional heater 
power K. The other parameter values of the system are c1=0.1, 
c=0.06, x;=0.3 and t= 0.2. 
3.2 阻尼 系数 的 影响 

为 了 研究 系统 中 阻尼 变化 对 系统 响应 的 影响 ， 
阻尼 模型 中 的 一 个 阻尼 参数 (cl ) 是 可 以 变化 的 。 
在 实际 的 实验 条 件 下 ， 系 统 阻尼 的 大 小 是 通过 改变 
管道 的 末端 条 件 来 改变 。 根 据 非 线 性 理论 ， 增 加 阻 
尼 可 以 增强 系统 的 稳定 性 。 图 6 为 系统 行为 演化 随 
PEAR Cc) 的 变化 的 分 兮 图 谱 。 如 预期 的 一 
样 ， 增 加 阻尼 确 具 有 增强 系统 稳定 的 作用 ， 因 为 平 
衡 解 对 于 较 大 的 阻尼 系数 ， 取 任意 时 间 延 迟 z 都 是 
稳定 的 ， 降 低 阻 尼 则 使 系统 失去 稳定 性 。 对 于 图 示 
的 时 间 延 迟 值 、 无 量 纲 加 热 功率 以 及 加 热 器 位 置 参 
数 ， 存 在 一 个 临界 值 = 0.19 ， 当 低 于 该 值 时 ， 系 
统 对 于 任何 幅度 的 初始 扰动 最 终 都 会 发 展 成 稳定 振 
幅 的 极限 环 ， 该 临界 值 即 为 亚 临 界 Hopf 分 贫 点 。 
Se, <0.19 时 ， 对 于 任何 幅度 的 初始 扰动 ， 系 统 最 
终 都 会 发 展 成 具有 稳定 振幅 的 极限 环 ; 当 
0.19<cl <0.25 时， 该 区 域 为 双 稳 态 区 域 ， 在 该 区 
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域内 ， 对 于 较 大 幅度 的 初始 扰动 ， 系 统 发 展 成 为 稳 亚 临 界 Hopf 分 岔 点 位 置 zm =0.14; 当 

定 极限 环 ， 而 对 于 小 振幅 的 扰动 ， 系 统 最 终 会 趋 于 0.14< x; <0.365 时 ， 系 统 处 于 全 局 不 稳定 区 域 ， 对 

平衡 ， 随 着 阻尼 系数 c, 的 进一步 增 大 ， 当 cl > 0.25 于 任何 的 初始 条 件 ， 系 统 将 渐 近 达到 相应 稳定 振幅 

时 ， 对 于 任意 幅度 的 初始 扰动 ， 系 统 都 最 终 趋 于 稳 的 极限 环 ， 当 0.365 < xy <0.39 时 ， 系 统 又 进入 双 稳 

定 状态 ，c = 0.25 即 为 折叠 点 。 态 区 域 ， 则 第 二 个 亚 临 界 Hopf 分 岔 点 位 置 
Xp_ = 0.365 ; 随 着 加 热 器 位 置 进 一 步 变 化 ， 当 


rs Xt > 0.39 时 ， 系 统 重 新 进入 全 局 稳定 性 区 域 ， 对 于 
1.6 4 : pe 
el 任何 幅度 的 初始 扰动 ， 系 统 最 终 都 会 趋 于 稳定 。 
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图 6 阻尼 系数 (c1) 变化 时 的 分 兮 图 ， 系 统 的 其 它 参 数值 —— “A, 一 一 一 一 “a -一 一 
X: KE1，c=0.06，xr=0.3，r= 0.2。 K 
Fig.6 Bifurcation plot for variation of damping coefficient c1. 
7 gy BE 变化 时 的 分 分 区 4 其 它 参 类 
The other parameter values of the system are K=1, c2=0.06, xy 7 加 热 器 位 置 Cp) RAR, RN ESS 
=0.3and r= 0.2. 为 : cl 0.1， C2 0.06, K=0.8, 7=0.2. 
3.3 加 热 器 位 置 的 影响 Fig.7 Bifurcation plot for variation of location of heater (x;). The 
加 热 器 的 位 置 (xp) 对 系统 动力 学 也 有 显著 的 other parameter values of the system are cl = 0.1, c2 = 0.06, K 


影响 ， 加 热 器 位 置 的 改变 是 通过 把 加 热 器 放置 在 沿 OS and r= 0.2. 
管道 长 度 不 同 的 位 置 处 来 实现 。 系 统 稳定 性 的 变化 ”3.4 时 间 延 迟 的 影响 
以 一 种 特别 的 方式 取决 于 管道 内 加 热 器 的 位 置 变 在 Rijke 管 中 ， 由 于 加 热 器 与 流动 介质 之 间 热 
化 。 加 热 器 位 置 Cy) 变化 时 的 分 岔 图 谱 如 图 7 所 IRTE (thermal inertia) 的 存在 ， 传 热 和 流动 速率 的 
示 。 当 管道 内 加 热 器 的 位 置 从 上 游 的 开端 开始 变化 变化 之 间 存 在 一 个 时 间 延 迟 r ， 其 大 小 与 加 热 器 中 
O ” 时， 该 系统 最 初 是 线性 稳定 的 。 在 加 热 器 的 位 置 xy ”电热 丝 的 直径 以 及 来 流速 度 有 关 。 在 实际 实验 条 件 
的 临界 值 ， 系 统 演变 为 线性 不 稳定 。 随 着 加 热 器 位 。 ”下 ， 时 间 延 迟 项 (+ ) 很 难 精确 的 进行 改变 。 时 间 
置 的 进一步 变化 ， 系 统 仍然 是 线性 不 稳定 的 ， 直 到 ER Cr) 变化 时 对 应 的 分 岔 图 谱 见 图 8。 亚 临界 
加 热 器 位 置 位 于 xzp， 在 该 点 形成 另 一 个 Hopf 分 岔 ”Hopf 分 兮 点 对 应 的 时 间 延 迟 值 为 r= 0.155 。 当 
点 ， 系 统 也 重新 变 成 稳定 状态 。 加 热 器 位 置 Q) z<0.13 时 ， 系 统 为 处 于 线性 稳定 状态 ， 对 于 任意 
ALA eZ, CERI, Hoof 分 命 点 均 为 ” ”幅度 的 初始 扰动 ， 系 统 最 终 都 趋 于 稳定 状态 ， 当 
亚 临界 Hopf 分 贫 。 在 两 个 Hopf 分 贫 点 之 间 ， 存 在 0.13<7<0.155 时 ， 系 统 对 于 小 振幅 的 扰动 最 终 趋 
一 个 全 局 不 稳定 区 域 。 在 全 局 不 稳定 的 区 域内 ， 对 于 平衡 ， 而 对 于 大 振幅 的 扰动 ， 系 统 发 展 成 为 极限 
于 任何 的 初始 条 件 ， 系 统 将 渐 近 达到 相应 稳定 振幅 。 ” 环 ， 即 双 稳 态 区 域 ， 当 +>0.155 时 ， 系 统 处 于 线性 
的 极限 环 。 根 据 加 热 器 位 置 变化 时 的 分 岔 图 ， 我 们 ” ”不 稳定 区 域 ， 对 于 任意 小 的 初始 扰动 ， 系 统 最 终 都 
可 以 得 到 两 个 亚 临界 Hopf 分 岔 点 的 具体 值 。 对 于 会 发 展 成 为 极限 环 。 值 得 注意 的 是 ， 发 现在 
我 们 选 定 的 其 他 系统 参数 ， 当 x < 0.11 时 ， 对 于 任 z> 0.155 的 线性 不 稳定 区 域内 ， 随 着 时 间 延 迟 〈= 
何 幅 度 的 初始 扰动 ， 系 统 最 终 都 会 趋 于 稳定 ， 当 ) DEREK, IRIE U, 的 值 先 逐 渐 增 大 后 又 逐渐 
0.11<xr <0.14 时 ， 系 统 处 于 双 稳 态 区 域 ， 系 统 根 。” ” 减 小 。 这 一 现象 可 通过 瑞 利 准则 来 解释 07。 时 间 延 
据 初 始 条 件 的 不 同 产生 不 同 的 演化 行为 ， 则 第 一 个 ”” 述 的 存在 是 因为 热源 的 放 热 和 流动 之 间 存 在 热 惯 


zy 


性 ， 当 加 热 器 功率 增 大 时 ， 时 间 延 迟 会 随 之 增 大 。 


时 间 延 迟 的 增 大 使 得 声场 振荡 和 放 热 脉动 的 相位 差 
发 生变 化 ， 当 二 者 相位 差 逐渐 增 大 而 远离 0” 时， 
振荡 幅度 会 逐渐 减 小 ， 从 而 形成 了 这 一 现象 。 据 此 


也 可 以 推测 当 r= 0.4 


时 ， 声 场 振荡 和 放 热 脉动 的 


相位 差 刚好 为 0”， 故 此 时 振荡 幅 值 最 大 。 


区 
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8 时 间 延 人 运 (r) 变化 的 分 岔 图 ， 系 统 的 其 它 参 数 为 


c=0.1, cs=0.06, K=0.8, xy=0.3。 
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4 结论 
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